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Abstract 

Growing interest is being given to transport of ultra-cold atomic gases through 
optical lattices generated by the interference of laser beams. In this connection we 
evaluate the phase-coherent transport of a spin-polarized gas of fermionic atoms 
along linear structures made from potential wells set in four alternative types of 
sequence. These are periodic chains of either identical wells or pairs of different 
wells, and chains of pairs of wells arranged in either a Fibonacci quasi-periodic 
sequence or a random sequence. The transmission coefficient of fermionic matter 
is evaluated in a T-matrix scattering approach by describing each array through a 
tight-binding Hamiltonian and by reducing it to an effective dimer by means of a 
decimation/renormalization method. The results are discussed in comparison with 
those pertaining to transport by Fermi-surface electrons coupled to an outgoing lead 
and by an atomic Bose-Einstein condensate. Main attention is given to (i) Bloch 
oscillations and their mapping into alternating-current flow through a Josephson 
junction; (ii) interference patterns that arise on period doubling and their analogy 
with beam splitting in optical inter ferometry; (hi) localization by quasi-periodic 
disorder inside a Fibonacci-ordered structure of double wells; and (iv) Anderson 
localization in a random structure of double wells. 
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1 Introduction 



Electronic transport in linear arrays of scatterers under an applied electric 
field has been extensively studied on grounds of both intrinsic fundamental 
interest and relevance to applications. A basic viewpoint was put forward 
by Landauer [1], who showed that the dc conductance is proportional to the 
phase-coherent transmittivity. From this viewpoint he proceeded to demon- 
strate that the resistance of a disordered one-dimensional (ID) lattice grows 
exponentially with the number of scatterers, as a result of wave-function lo- 
calization. Within a tight-binding framework the scatterers in the array can 
in fact be taken to vary according to a chosen law (such as periodic, quasi- 
periodic, or random) and the density of states and the transmission coefficient 
can be calculated by Green's function methods (see for instance Economu [2]). 
A simplifying assumption treats each scatterer as a potential well, thus omit- 
ting a detailed account of its internal structure but allowing analytical results 
to be obtained for linear arrays (see for instance [3,4,5]). 

Advances in the preparation and manipulation of ultra-cold atomic gases have 
very significantly broadened the scope of the study of phase-coherent quan- 
tum transport through arrays of potential wells (for a recent review see [6]). 
A standing-wave created by the interference of two counter-propagating laser 
beams, which are detuned away for an atomic absorption line, can be su- 
perposed onto an elongated magnetic trap to realise an almost ideal linear 
periodic array of potential wells for an atomic gas. Higher-dimensional arrays 
have been created by using two or three pairs of lasers. Experiments on Bosc 
gases in such "optical lattices" have led to the observation of Bloch oscilla- 
tions under the force of gravity both in ultra-cold atomic gases [7] and in 
a Bose-Einstein condensate [8], Landauer-Zener tunneling [9], Josephson-like 
oscillations [10,11], a superfluid to Mott-insulator transition [12], and the ID 
band structure [13]. Collective Bloch oscillations lasting for very many peri- 
ods have been observed in a spin-polarized Fermi gas inside an optical lattice, 
opening a novel route to high-precision interferometry and to the measure- 
ment of forces with microscopic spatial resolution [14]. The ideal conditions of 
extremely low temperature and high purity in which experiments on atomic 
gases can be carried out offers a unique opportunity for the investigation and 
testing of transport theory at unprecedented levels of accuracy and depth. 

A Bose-Einstein condensate that has been adiabatically loaded into an array 
of potential wells keeps its long-range phase coherence, provided that the bar- 
riers between neighboring wells are not too high and still allow tunnel between 
them [15,16]. In this situation a weak external force accelerates the conden- 
sate through the band states of the array as if it were a single quasi-particle 
[17]. In previous work we have treated coherent transport of an atomic con- 
densate through a variety of linear arrays, using a decimation/renormalization 
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method in a tight-binding approach to reduce each array to an effective dimer 
[18,19,20]. The arrays that we considered were of four basic types, i.e. periodic 
chains of either identical wells or pairs of different wells and chains of pairs 
of wells arranged in either a Fibonacci quasi-periodic sequence or a random 
sequence. In each case we also proposed a specific set-up of optical lasers that 
through beam interference could create the array in the laboratory. 



In this work we extend the theory to a spin-polarized Fermi gas, this expres- 
sion being a standard short-hand notation to indicate that the fermionic atoms 
occupy a single Zeeman sub-level in a magnetic trap. Spin-polarized fermions 
in an ultra-cold gas can to a very good approximation be regarded as non- 
interacting, since the antisymmetry of the many-body wave function under 
exchange that is embodied in the Pauli principle inhibits close encounters be- 
tween them and this suppresses s-wave collisions (for a full discussion see for 
instance [6]). The absence of interactions makes the spin-polarized Fermi gas 
a good candidate to observe Anderson localization of matter waves. In the 
ground state the fermions occupy the single-particle levels of the confining 
potential up to a maximum fixed by their number, and in particular in a ID 
single-well periodic lattice at half filling they fill the bottom half of the lowest 
energy band. This situation should be contrasted with that of a Bose-Einstein 
condensate, where the gas in its ground state is condensed in the lowest (zero 
quasi-momentum) band state and, as already noted, is driven by a constant 
external force through the band states as if it were a single quasi-particle. 
Coherent motion of a spin-polarized Fermi gas under an imposed bias can 
still occur if each fermion is in turn accelerated into successively higher band 
states. The whole distribution is thereby coherently accelerated through the 
band structure in the absence of dissipative scattering. 



The contents of the paper are briefly as follows. In Section 2 we present the 
model and discuss the essential aspects of the density of states and of the 
calculation of the transmission coefficient. Numerical results are reported in 
Section 3 for a spin-polarized gas of fermionic 40 K atoms, this being the system 
for which long-lasting Bloch oscillations in a single-well optical lattice have 
been observed by Roati et al. [14]. Our results are discussed in comparison 
with those for coherent transport by Fermi-level electrons [5] and by a Bose- 
Einstein condensate [18,19,20]. The main points at issue are (i) the Bloch 
oscillations of the gas and their mapping into current flow across a Josephson 
junction; (ii) the interference patterns between Bragg scattered inter-subband 
tunneling fragments that arise on period doubling; (iii) the localization of 
fermionic matter induced by the introduction of quasi-periodicity; and (iv) 
Anderson localization induced by randomness. Finally, Section 4 summarizes 
our conclusions. 
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2 The tight binding model and the array transmittivity 



We consider a ID sequence of iV equally spaced potential wells to be occupied 
by Nf spin-polarized fermions. In a tight binding model the Hamiltonian is 

N 

H = E( e * I *) <* I + (ki+i + *<+i,il * + !) (» I )} W 



where the site energy corresponds to the lowest energy level in the i-th. 
well and t i}i ±i are the hopping energies between adjacent wells. The case of 
practical interest for an atomic gas in a ID optical lattice is that in which a 
few handred sites are occupied. The calculation reported in Section 3 below 
will refer to 40 K atoms at half filling, with N f = 100 and N = 200. 

The array is described by an external potential U(z) constructed from the 
periodic function sin 2 (irz/d), d being the distance between adjacent wells, and 
a bimodal distribution of well depths (Ui and U 2 , say). The potential created 
on an atomic gas by a standing optical wave generated from two counter- 
propagating laser beams correspond to JJi = U\ and d equal to one-half the 
laser-light wavelength. In addition to this case, in which U(z) describes a 
period sequence of identical wells, we also treat cases in which the C/j's are 
chosen so that the well depth (i) regularly alternate along the array, or (ii) 
form a Fibonacci sequence, or (iii) form a sequence produced by a random- 
number generator. 

The global density of states of each array is obtained from the Green's func- 
tion of the model by means of the Kirkman-Pendry relation [21], expressing 
coherence between the first and the last site of the chain. Illustrative examples 
have been given in our previous work [5]. Of course, period doubling breaks the 
lowest energy band of the chain into two sub-bands separated by an energy 
gap, whereas the introduction of quasi-periodic or random disorder induces 
spectral fragmentation and a multiplicity of pseudo-gaps. 

An applied constant force causes a linear tilt of the spectrum and the atoms 
travel along the chain and at the same time are accelerated through the energy 
spectrum [22]. A schematic representation of the trajectory of an atom moving 
in space and energy can be obtained from the local (site projected) density of 
states Di(E) (LDOS) and is shown by the arrows in Fig. 1 for a chosen value 
of the initial energy. The LDOS is given by the expression 

D i (E) = --Im(i\(E-H)- 1 \i) =--ImG iti (E) (2) 
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Fig. 1. Local density of states as a function of position and energy for a single- well, 
a periodic double-well, a Fibonacci, and a random array (from top left to bottom 
right). The arrow indicates the trajectory of an atom starting at a given initial 
energy. 



where Gi^(E) is the site-projected Green's function at energy E+ir] in the limit 
of vanishing positive r\. This matrix element can be evaluated as G^i(E) = 
[E — ei(E)]" 1 by reducing the array to an effective site with energy Si(E) 
through decimation of all j-th sites with j ^ i, the renormalized Hamiltonian 
being a c-number expressed in term of a continued fraction [23,24,25,26,27]. 
The LDOS in Fig. 1 shows that at a given energy the population of a subset of 
sites is favored in both the double-well and the disordered arrays, these sites 
being rather irregularly distributed along the array in the latter case [19]. 

The same renormalization procedure is used in the calculation of the transmit- 
tivity to reduce each array to an effective dimer containing just the first and 
the last site with renormalized site and hopping energies [23,24,25,26,27]. The 
Hamiltonian H(E) of the dimer depends on the energy E and is expressed as 
a 2 x 2 matrix, 



H(E) 



(3) 



Here, the effective site and hopping energies are determined by the recursive 
relations 
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e^(E)=er i> (E) + t w (E) - . x t, +1 , +2 , (4) 

/ — r v,; 

1 



S-+2(^) - e i+2 + <j+2j+i — — iy— , (5) 



(#) - iij+i (£■) — 1(— iy— * j+i.j+2 , (6) 



and tj + i t i = t±j + i for j > 1. The initial values are given by the original 
Hamiltonian parameters, e^ \E) = e$ and tip(E) = tip. 



2.1 The transmittivity 



We can now proceed to calculate the transmittivity of the array for a gas 
driven by a constant force F. For this purpose we calculate the scattering 
matrix in the presence of the bias, in a configuration in which the array is 
connected to incoming and outgoing leads [18]. The leads mimic a coupling to 
the continuum by injecting and extracting a steady-state particle current. In 
describing an experiment such as that carried out by Anderson and Kasevich 
[8] on a Bose-Einstein condensate, the incoming lead would correspond to a 
continuous replenishing of the condensate in the array and the outgoing lead 
to a detecting system counting the atoms that leave the array in coherent 
drops separated in time by the Bloch-oscillation period under gravity. In an 
experiment such as that carried out by Roati et al. [14] on a Fermi gas, instead, 
the leads should be blocked after loading the gas in the array. 

With specific reference to a spin-polarized atomic Fermi gas, all fermions with 
initial momentum p m in the range p- m G [— Pf,Pf\, with pp = (2?71-Ef) 1//2 being 
the Fermi momentum, move coherently under the applied force in the absence 
of dissipative scattering processes. We define iV ou t = (E m3iX — E in ) / (Fd) (with 
^Vout < N) as the number of hops after which an atom having initial energy E in 
and momentum p in reaches the highest point in the dispersion curve at energy 
-E'max and can either leave the array through the outgoing lead or continue its 
trajectory through the band states (see Fig. 2 for an illustration referring to 
the periodic single- well and double- well arrays). The leads are defined by two 
additional terms in the Hamiltonian, 

#L,in= E {4 n |n)(n|+t (|n)(n + l|+c.c.)} (7) 

n=— oo 
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p=h/2d 




Nout 




Nout 



Fig. 2. Band tilting in space for a single-well array (left) and a periodic double-well 
array (right), due to a constant force F. Coupling to an incoming and an outgoing 
lead is also shown for an atom with initial momentum p- m and exiting after N ont 
hops at the Brillouin zone edge (left) or at the Brillouin zone center (right). The 
position z is in units of d. 

and 



#L,out — 



oo 

£ 

n=N+l 



{E° ut \n)(n\ +t (\n- l)(n| +c.c.)}. 



We take the hopping energy t Q equal to one fourth of the spectral band width 
and shift the centers of the leads by FdN out in order to optimize their coupling 
to the array. 

The outgoing wave function |</J„ ut ) is an eigenstate of the outgoing lead and 
is related to the incoming wave function |y2f n ), which is an eigenstate of the 
incoming lead, by 

lv4t> = (i + G°r)lvO- (9) 



The wave vectors k and k of the two wave functions are uniquely determined 
by the relations 

1 fE-E^ out \ , . 

k, k = — arccos . (10) 

d \ 2t n 



Further, in Eq. (9) G° = (E — Hq)^ 1 is the Green's function of the disconnected 
leads, described by the Hamiltonian H = H L in + H L out , and the scattering 
matrix T(E) = H\(l — G^Hi)^ 1 referred to the perturbation Hamiltonian 

Hj = H\ N=Nout - (4 n + £o° ut l^out)(iV out |), (11) 

with H\N=N ou t being given by Eq. (3) for the case N = N out . 

Since the wave functions \(p? n ) and |<^o Ut ) are denned in disconnected spaces, 
the projection of |y?out) onto the localized function |n) can be written as 
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(n|¥4t>=E^,i r wO'K> 

3<l,l 

= V2G% ou T NontA sm(Kd)u K (d) (12) 

where we have set (1 \ip^ n ) = u K (d)(e lKd — e~ lKd ) / (i\/2) with u K (d) being the 
Wannier function in the potential U(z). The Green's function element G° iA r out 
in Eq.(12) determines the coherence between site n and site N out on the chain 
for the outgoing lead and can be written as 

rn-N out 

G°n, Nout = ..n+l-^i eX P[^ «(« + 1 - ^ut) d] . (13) 



In the out-of-equilibrium picture, the velocities fj n and v out of the incoming 
and outgoing wave functions enter the definition of the transmission coefficient 
T as 



lim n ^ +00 (n! 


¥>OUt)(^OUt 


I^Kut 


lim^.ooH^fJ^fJ 


m)v- m 



(14) 



Finally we obtain 



T = 4 ' ri '^ out 1 sin («d) sin (Rd) . (15) 
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We remark that Eq. (15) reduces to the transmission coefficient evaluated in 
[5] for Fermi-level electrons when we set Eq u1 = E n . 



3 Numerical results 



Our task now is to determine the parameters and tj^+i entering the Hamil- 
tonian in Eq. (1). We proceed for this purpose to a ID reducing of the system 
by introducing the transverse width a± of the atomic cloud in an elongated 
cigar-shaped harmonic trap and the ID wave function in the i-th well. In a 
tight binding scheme the Wannier function ^i(z) in the potential U(z) can be 
written, according to the work of Slater [28], as a Gaussian function having a 
width a Zi determined by the harmonic approximation to the well. That is, 

= (ttVVJ- 1 exp[-(* - z i Y/(2al)} (16) 
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where a Zi = (h/muji) 1 ^ 2 with uji = (4f/ i -E R ) 1 / 2 , E R being the recoil energy. The 
site energies are then given by 



2V72 



2m 



+ U(z) + maz + C 



(17) 



Here, a = F/m is the acceleration due to the constant force F applied to the 
fermions, and C = h 2 j (2ma 2 _) + \muj\_a\_ is a quantity taking into account 
the reduction to ID dimensionality with u± being the radial frequency of the 
trap [29]. The hopping energy are similarly given by 



U,i+i = J dz^i(z) 



2m 



+ U(z) + C 



(18) 



For a numerical illustration we have chosen to adopt typical system parameters 
used in experiments on atomic gases, by setting U\ = 3.5Er and 2d = 763 nm 
with (U 2 - Ui)/U! = 10 -2 in the double-well arrays. Our results are shown in 
Fig. 3 where we have plotted the transmittivity of the spin-polarized Fermi 
gas of 40 K atoms through periodic arrays and through a Fibonacci-ordered 
and a random array of double wells, as a function of the inverse of the applied 
force and of the energy. In all cases the array is at half filling and the inverse 
force is described through the quantity Ts/n, where Tb = h/{2Fd) is the 
Bloch-oscillation period for the fermions in the double-well lattice and r t is a 
time-scale parameter to be defined immediately below. The energy variable is 
written as the ratio (E — E min )/(E F — E min ). 

It is seen from Fig. 3 that fermions at the Fermi level in the single-well lattice 
have a better coupling with the outgoing lead at low values of the drive, while 
those at very low energies have a greater coefficient of transmission at large 
values of the drive. The behavior of the transmittivity through the double- 
well lattice is illustrated in the top-right panel of Fig. 3. We expect in this 
case destructive interference in the transmission whenever the Bloch period 
is an integer multiple of the time to tunnel twice across the inter-subband 
minigap AE. The underlying phenomenon is the interference occurring at 
p = —h/(Ad) (the left-hand boundary of the Brillouin zone in the split-band 
configuration) between the atomic wave packets that are Bragg scattered at 
p — h/ (4d) in the lower sub-band and those that have tunneled into the upper 
sub-band, traveled through it, and tunneled back into the lower sub-band at 
p = —h/(4d). From our numerical data we determine the fermion tunneling 
time as r t = 6irh/AE. A proportionality of the tunneling time to h/AE is 
expected from basic quantum mechanics (see for instance [30]), but the specific 
value that we find for the proportionality constant naturally depends on the 
model Hamiltonian that we have assumed. 
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T 



0.1 





Fig. 3. Transmittivity T (in log scale) of a spin-polarized Fermi gas through 
a single- well lattice (top-left), a periodic double- well array (top-right), a Fi- 
bonacci-ordered double- well array (bottom- left), and a random double- well array 
(bottom-right). The transmittivity at half filling is plotted as a function of Ts/rt 
and of (E — E ni i n )/(Ep — E m [ n ). The time scales Tb and tj in all cases are taken as 
the Bloch oscillation period and the inter-subband tunneling time in the double-well 
lattice. 

Turning to the effect of quasi-periodic or random disorder, we first remark 
that the fermions traveling through the array of Fibonacci-ordered double 
wells explore an energy spectrum which on average is rather more akin to 
the single-well band structure than to the double-well one, but is regularly 
fragmented by a multiplicity of pseudo-gaps [5]. The quasi-periodic disorder 
nevertheless induces marked peaks and troughs in the transmission as can 
be seen in the bottom-left panel of Fig. 3. The peak positions do not show 
any regularity and depend on both the applied drive and the fermion energy. 
Finally, the bottom-right panel in Fig. 3 shows that random disorder in a 
double-well array at half occupancy induces an interference pattern which is 
qualitatively similar to that of a Fibonacci array. 

With the aim of making a direct comparison of our results for transport in a 
Fermi gas with the transmittivity of a Bose-Einstein condensate calculated in 
previous work [18,19,20], we have plotted in Fig. 4 the mean transmittivity of 
the Fermi gas after averaging over incoming states, namely the quantity 




f =[ J T(E)D in (E)dE}/[ J D in (E)dE] 



(19) 
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Fig. 4. Averaged transmittivity f of a spin-polarized Fermi gas at half filling as 
a function of Tb/t^. The various panels refer to a single- well lattice (top- left), a 
periodic double-well array (top-right), a Fibonacci-ordered double- well array (bot- 
tom-left), and a random double- well array (bottom-right). 



as a function of Tb/tj, D- m (E) being the density of states of the incoming lead. 
The behavior of the averaged fermion transmittivity is qualitatively similar to 
that found for bosons. In particular, Bloch oscillations with a period h/(Fd) 
are executed by the whole distribution of fermions once the coupling to the 
leads is cut. This type of behavior was observed for a condensate in the ex- 
periments of Anderson and Kasevich [8] and was interpreted in term of the 
condensate moving through the band states as a single coherent quasi-particle 
[15,18]. Our calculations confirm that a similar coherent-transport behavior, 
as observed in the experiments of Roati et al. [14], is exhibited by a spin- 
polarized Fermi gas in a single-well lattice, even though the fermions start 
from a ground-state configuration in which they are distributed over many 
band states up to the Fermi level instead of being all condensed in the lowest 
band state at the Brillouin zone center. 

However, it is also seen from Fig. 4 in comparison with previous results for 
a Bose-Einstein condensate [18,19,20] that inter-subband coherence is sharply 
reduced in a Fermi gas and that the consequences of quasi-periodic or random 
disorder in a double- well array are more dramatic. The low values in the 
transmission coefficient in the presence of disorder essentially correspond to 
Anderson localization of the fermionic matter induced by the large number of 
pseudo-gaps that have opened up in the density of states. The non-zero values 
of the transmittivity are due to finite-size effects. 
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4 Conclusions and future perspectives 

In summary, we have focused attention in this work on some main features 
of the coherent transport of gaseous ultra-cold fermionic matter through lin- 
ear arrays of potential wells under a constant external drive. Our present 
results reduce in the case of infinitesimal bias to those holding for transport 
by Fermi-level electrons through arrays of confining sites such as quantum 
dots, as evaluated in previous work [3,4,5]. The discussion of our results has 
been framed within a comparison with those for coherent transport of bosonic 
matter by a Bose-Einstein condensate [15,16,17,18,19,20] with main attention 
to long-lasting Bloch oscillations in a lattice of confining sites, to matter-wave 
interference induced by splitting through an inter-subband minigap and hav- 
ing an analogue in optical beam-splitting experiments, and to localization of 
matter waves by quasi-periodicity and by random disorder. We have used a 
procedure of site decimation followed by renormalization, which allows ana- 
lytic development of the theory for linear arrays and could usefully be extended 
to transport in networks of confining sites having higher dimensionality. 

In comparison with the many experimental studies that have been carried out 
on electronic transport through arrays or networks of quantum dots, the ex- 
perimental study of matter transport by quantum atomic or molecular gases 
through optical arrays is still in its infancy though promising of novel funda- 
mental developments. Both bosons and fermions have been observed to per- 
form Bloch oscillations inside a single- well optical lattice [8,14] and in both 
cases the experiments give clear evidence for long-range coherence of matter 
waves. Transport of quantum matter through other types of array remains to 
be studied experimentally. A most promising development would seem to lie 
in the direction of one-dimensional transport by two-component Fermi gases 
consisting of two spin populations, forming a Luttinger liquid for repulsive 
interspecies interactions and leading to Luther-Emery spin pairing for attrac- 
tive interactions [31]. Separate control of the two spin populations has been 
achieved [32,33] and opens the way to the study of exotic superfluid states 
[34,35,36,37,38]. 
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